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Abstract. Twist stiffness and an asymmetric bending stiffness of a polymer or a polymer bundle is captured 
by the elastic ribbon model. We investigate the effects a ring geometry induces to a thermally fluctuating 
ribbon, finding bend-bend coupling in addition to twist-bend coupling. Furthermore, due to the geometric 
constraint the polymer's effective bending stiffness increases. A new parameter for experimental investiga- 
tions of polymer bundles is proposed: the mean square diameter of a ribbonlike ring, which is determined 
analytically in the semiflexible limit. Monte Carlo simulations are performed which affirm the model's 
prediction up to high flexibility. 
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1 Introduction 



j— ( \ The semiflexible polymers constituting the cell's cytoskele- 
t-H ■ ton determine morphology and elasticity of the cell. By use 
' of binding proteins the polymers organize into bundles of 
00 . closely packed filaments, which form the building blocks 
' for a variety of cellular processes such as filopodia, villi 
[--■« or cilia [1.2,3,4 . Recent investigations of the mechanical 
■ properties of polymer bundles [5] pose the question for a 
| coarse-grained description of their mechanics beyond the 
• <-h . macroscopic models for semiflexible polymers such as the 
' Kratky-Porod model [5] and its continuous description [7] , 
^ \ which represent the polymer by a space curve with a sin- 
- - ' gle material parameter. A more general description is ob- 
tained by the classical elastic ribbon [5], which accounts 
for three parameters concerning bending and twisting. 
There are also advances to describe polymer bundles as 
a set of interconnected semiflexible polymers, starting out 
with railway track models of two coupled filaments in two 
|9J and three [10 l fTl ] fT2j dimensions, whose exact imple- 
mentation and results were much in debate [12) . Recently, 
a model for wormlikc bundles with isotropic cross-section, 
consisting of a finite number of semiflexible polymers, has 
been introduced and analyzed [13lfT4] . These models show 
new material properties, such as length-dependent bend- 
ing and twisting moduli [13], as they include stretching 
and shearing terms, not accounted for in the elastic ribbon 
model considered in this paper. The latter regards both 
bending and twist moduli as fixed material parameters. 

Nature not only imposes geometric constraints on bio- 
polymers by confinement through membranes but also de- 
liberately uses the advantages of certain geometries as in 
the case of circular DNA. Especially the ring geometry 



induces interesting effects. In the context of plectonemes 
mechanical equilibria of elastic rings have been studied 
extensively [15j . whereas thermal equilibria have mainly 
been addressed for flexible rings pT6 | [TT | fT8] . Only recently 
investigations dealt with the thermal fluctuations of poly- 
mer ribbons that are stress-free in a circular configuration 
19J. However, both the mechanical and the thermal mo- 
tion of rings with intrinsic bend differ principally from ini- 
tially straight polymers as has been observed for the plec- 
tonem transition [20 . For biopolymers participating in 
the cytoskelcton such as actin and microtubuli the stress- 
free conformation is straight. Even small DNA rings have 
lately been found to exhibit no intrinsic bending without 
attachment of proteins [21 . Furthermore, ribbonlikc rings 
require theoretical modeling as in vitro experiments have 
shown that the mechanical properties of polymer bundles 
can be well studied when the bundle is constrained to a 
ring structure 5 . In vivo polymer bundle rings are found 
in erythrocytes of birds and reptiles [22] . 

Ligating a ribbon's end to form a ring reduces the 
set of conformations available to the ribbon. Perturbing 
the ring's mechanical equilibrium configuration by a small 
twist induces a bending of the ribbon's center line, as 
can be easily visualized with a ring consisting of a small 
strip of paper. This twist-bend coupling is mathematically 
captured by White's formula 23 , 24 , which connects the 
overall twist to the global configurational integral denoted 
writhe for any closed curve. However, we will show that 
for a ribbonlike ring with asymmetric cross section bend- 
ing within the equilibrium plane and bending transversal 
to the equilibrium plane are also coupled. In this case it is 
instructive to parameterize the ring in terms of Euler an- 
gles such that each of them describes either a bending or 
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a twisting motion. Modeling thermally fluctuating rings 
we pursue to give analytical expressions for experimen- 
tally observable properties of the ring. In analogy to the 
mean square end-to-end distance of a polymer the mean 
square diameter is well suited to represent the statistical 
properties of a fluctuating ribbonlike ring. 

In this work we investigate the effects of a ring geome- 
try on a fluctuating semiflexible polymer ribbon. In Sec- 
tion [2] the elastic ribbon model is formulated in the limit 
of small fluctuations about a ring. The elastic free energy 
already indicates the first effect of the ring geometry on 
a ribbon: the bend-bend coupling in addition to twist- 
bend coupling. Furthermore, in Section^ the foundations 
are laid for the analytic calculation of the experimentally 
accessible mean square diameter. The discussion of this 
ensemble parameter in Section [3] reveals the quality of the 
coupling and the second effect due to the geometric con- 
straint: the effective stiffening of the ribbon. Assessing the 
quality of the analytical result for the mean square diame- 
ter with Monte Carlo simulation we find that our model 
predicts the correct behavior for a polymer with symmet- 
ric cross section up to total contour length seven times the 
persistence length. Beyond this stiff limit our simulations 
confirm previous results for almost flexible polymer rings. 
Thus, polymer rings with symmetric cross section are fully 
characterized analytically. 



2 The model 

To specify a deformed state of a ribbon a local body co- 
ordinate system {ti(s), t 2 (s), t 3 (s)} is assigned to every 
point s € [0, L] along the center line of the ribbon, such 
that the t 3 (s) axis points tangent to the center line of 
the rod in the direction of increasing s, while ti(s) and 
t 2 (s) are aligned to the principal axes of the cross section. 
The change of body coordinates along the inextensible rib- 
bon is described by rotations ui via the generalized Frenet 
equations: 

^ =u xU, i = 1,2,3. (1) 
as 

u)\ and ll>2 may be interpreted as the curvatures in the 
principal directions, while cl> 3 may be denoted as the heli- 
cal deformation density [15]. Hence, the elastic free energy 
of a ribbon reads [8j: 

F = —— I ds [axujj + a 2 io\ + a 3 w|] , (2) 
z Jo 

where ai and a 2 denote the bending stiffnesses along the 
principal axes of the cross section and et 3 indicates the 
twist stiffness. For the geometry of a ring it is instruc- 
tive to express the rotations uj in terms of Euler angles 
4>(s), 0(s), ip(s) according to Euler's equations, where cp is 
the azimuthal angle, 9 the polar angle and ip the twist 
angle. The constraint of a closed ring of contour length 
L = 2irR c is then satisfied by 4>o(s) = s/R c and &o( s ) = 
7r/2. The remaining angle ipo{s) is chosen as to minimize 
the elastic free energy which is achieved by ipo(s) = n/2 




Fig. 1. Illustration of the Euler angles and their coupling in- 
duced by the ring geometry: A sinusoidal variation of each 
Euler angle at a time results in more than one free energy 
term, if coupling is present: (a) 5<p — sin(s/-R c ), in-plane bend- 
ing only, F = a%, (b) 88 — sin(s/i? c ), transversal bending 
yields in-plane bending and twisting, F = — ai + a,2 + 0,3, (c) 
dtp = sin(s/R c ), twisting yields both in-plane and transversal 
bending, F — —ai + 0,2 + 0,3; F in units of ttIcbT /2R c . 



for a\ < et 2 or by ipo{s) — for a\ > a 2 . These two cases 
are essentially the same, as the major principal axis of the 
cross section is always perpendicular to the equilibrium 
plane of the ring. We therefore restrict our discussion to 
ai < a 2 . As the biopolymers constituting polymer bundles 
are itself semiflexible, we investigate the limit of small fluc- 
tuations of the ribbonlike ring. In this limit self-avoidance 
is satisfied and hence our phantom chain model should 
yield correct results for real polymers. Furthermore, we 
keep the total length of the ribbon strictly fixed. Small 
fluctuations in the contour length might though lead to 
interesting behavior, the investigation of which lies be- 
yond the scope of this publication. 

Assuming small deviations Sep, SO, Sip from the ground 
state of the ring specified by (po,0o,ipo the rotations u> 
are expanded and all quadratic terms obtained from the 
squared rotations are taken into account to obtain the 
elastic free energy in the limit of small fluctuations: 

I-, k B T r L J \(d8cj>V fse\\fdse\ 2 ' 

F =—J ds \ ai [{-d7) -{rJ + UrJ_ 

. ( Sip dse\ 2 fdSip se\ 2 \ , Q , 

In parameterizing the ribbon in terms of Euler angles a 
distinct motion was assigned to each angle. The azimuthal 
angle Sep describes bending within the equilibrium plane, 
while the polar angle S9 characterizes transversal bend- 
ing and the twist angle Sip depicts twisting around the 
center line of the ribbon, see Fig. [TJ Considering the elas- 
tic free energy associated with each conformation of the 
ribbon, three energy components can be distinguished, the 
in-plane bending energy -Fin-piano, the transversal bending 
energy Ftrans and the twisting energy F tW i S t, recognizable 
by their proportionality to their associated stiffnesses a\ , 
a 2 and a 3 , respectively. In a completely decoupled sys- 
tem such as an open polymer ribbon each energy com- 
ponent would comprise just the spatial derivative of its 
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corresponding Euler angle 



b = I ds< a\ 



\ ds 



fdSipY 

fl3 UrJ 



(4) 



In the limit of R c — > oo, keeping a±, a 2 and 03 fixed, the 
elastic free energy of a polymer ribbon ring as given in 
equation decouples to the above expression, as local 
fluctuations are not influenced by the ring geometry any- 
more. However, for small contour radius with respect to 
the bending and twisting stiffnesses, the elastic free energy 
of a polymer ribbon ring exhibits both twist-bend coupling 
and bend-bend coupling. Twist-bend coupling heuristi- 
cally means that twisting induces bending and vice versa. 
In the case of F tw i s t the twist arising form the change of 
the twist angle Sip is diminished by the height 59 relative 
to the contour radius R c the ring gains through transver- 
sal bending. The opposite coupling effect is observed for 
-Ftrans where the bending energy is augmented by the lo- 
cal twist relative to R c . Analyzing the in-plane bending 
energy component 



in — plane 



(d5f\ 2 _ fse_ s 2 



\ ds J 



Rc 



dS6 6tfj 
J ~ds~~R r 



(5) 



twist-bend coupling is found again, but now also a varia- 
tion of the transversal bending Euler angle 59 reduces the 
in-plane bending energy. This is identified as bend-bend 
coupling. Furthermore, the release by the twist depends on 
the curvature of the transversal bending, hence, a mixed 
term arises. 

As a demonstrative example the conformations of a 
ribbon with sinusoidal variation of each Euler angle are 
depicted in Fig. [TJ Both altering 68 and Sip causes contri- 
butions to all three terms in the elastic free energy exem- 
plifying the couplings. Remarkably, altering the twist an- 
gle 6tf) or the polar angle 69 in the same sinusoidal manner 
requires exactly the same energy cost. Furthermore, for a 
ribbon with symmetric cross section the only energy con- 
tribution due to sinusoidal transversal bending or twisting 
is proportional to the twist stiffness only, as other terms 
cancel. Previous work assumed the unstressed state of a 
polymer to be a ring conformation |I9j . whose applica- 
bility to cytoskeletal polymers maybe limited, yielding no 
coupling of transversal bending and twisting to in-plane 
bending. 

Given the elastic free energy statistical properties may 
be derived. Because of its importance for experimental in- 
vestigations we shall focus on the mean square diameter of 
the ring. Since the Euler angles must obey periodic bound- 
ary conditions to satisfy a closed ring they are expanded in 
a Fourier series. In the Fourier representation the periodic 



boundary condition correspond to: 59(0) = S(fi(l) = 0. 
The Fourier transformed elastic free energy then results 



F = 



Rc 



{(oa-aOI^O)! 1 



+ (a 3 + a 2 - ai )\i59{l) + 5^(l)\ 2 

OO 

[("2 - a{)\m86{n) + 6ip(n)\ 2 + ai n 2 \5^{n)\ 2 
(n 2 - l)\69(n)\ 2 + a 3 |m#(n) - 59{n)\ 2 ] } . (6) 



n=2 



Note that concerning the Euler angles themselves only 
the twist angle Sip and the polar angle 69 are coupled. 
Diagonalizing the Fourier transformed elastic free energy 
and applying the equipartition theorem the correlations 
of the Euler angles are obtained depending on the mean 
square modes. 

7? °° 

(<ty(*a)ty(«i)> = — E(^ 2 W> cos(ns/i? c ) , (7) 

IT ' * 



(56(s 2 )66( Sl )) 



R c cos(s/R c ) 
it a 3 + a 2 — ai 

+— Y(69 2 (n)) cos(ns/R c ) , (8) 

7T * — ' 

R c , R c cos(s/R c ) 



7r(a 2 — ai) n a 3 + a 2 — a± 

7T 



^(<5V; 2 (n))cos(Wi? c ) , (9) 



where s = |s 2 — si|. Closed expressions for the modes 
are only obtained for n > 2. For the azimuthal angle 5<p 
the zeroth mode vanishes due to the rotational symmetry 
around the axis through the ring's center of mass. The 
first mode is zero due to the periodic boundary conditions, 
which also set the zeroth polar angle mode equal nought. 
However, the zeroth mode of the twist angle Sip vanishes 
only in the limit of (a 2 — ai) — » 00. The minimum of the 
elastic free energy is given by 6ip = 0, such that the major 
principal axis of the cross section along a 2 in our definition 
is perpendicular to the equilibrium plane. For finite val- 
ues of both bending stiffnesses small fluctuations around 
5ip = arise with an amplitude that increases with a de- 
caying difference between the bending stiffnesses, being 
not defined for equal magnitude. Among all three stiff- 
nesses a\ influences the fluctuations of a polymer ribbon 
ring the most. For a\ = 00 all deformations of the ring are 
energetically hindered, since also twisting and transversal 
bending require suppressed in-plane bending due to the 
coupling. For small a\ on the other hand the fluctuations 
exceed the small deviations approximation. As the effect 
of the other two parameter is smaller it is reasonable to 
study their significance relative to the in-plane bending 
stiffness a\. Hence, we introduce the dimensionless pa- 
rameters a = ai/a 2 and r = 01/03, in which the modes 
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of the Euler angles for n > 2 are given by: 



(S0 2 (n)) = 
WV)> = 



Rc 

01 
Oi 



1 



1 



- 1 



1 



+ (l-r)(a-l) 



+ (l-r)(a-l)J 



(10) 
(11) 
(12) 



As variations in 5<f> do not influence transversal bending or 
twisting, the modes only depend on the in-plane bending 
stiffness a%, whereas the modes of the multiply coupled 
Euler angles 69 and 6cj) are functions of all three stiffnesses. 

The mean square diameter is by definition just the 
mean square distance between two positions r(s) sepa- 
rated by half the contour length L/2 = ttR c along the 
center line of the ribbon. The mean square distance itself 
can be calculated from the tangent-tangent correlation by: 

CD 2 ) = <(r(L/2)-r(0)) 2 ) 



dy / dy'(t z ( V )hW)) ■ (13) 
'o Jo 

The tangent is approximated for small deviations from the 
rigid ring up to second order in 8(f> and 59 such that the 
tangent-tangent correlation yields: 



{h{si)h{s2)) = (Scf)(sx)6<f>(s2)) cos 
+ (69( Sl )S9(s 2 )) 
+ [l - (59 2 ) - (5<t> 2 )] cos 



|fl ~ 32 | 

R, 



Sl - S 2 



Rc 



(14) 



Evaluating the double integrals the mean square diameter 
is represented in terms of the modes of the Euler angles: 

AT? 2 °° AR 2 °° 

(D*) = {2R c f - ¥k J2(54> 2 (n)) - ^ X>~ 2 (n)) 

n=2 n=2 

TR| ~ (n 2 + l)(^») 
7T ^ (n 2 -l) 2 

n even 



4i? 2 



OO 

£ 

n=3 
n odd 



(69" (n)) 



(15) 



From this expression the mean square diameter can be de- 
rived. The sum of all four sums always yields a negative 
term corresponding to a decrease of the mean squared di- 
ameter from its rigid ring value of D = 2R C . Note that 
the mean square diameter does not depend on the twist 
modes since twist motions in the weakly fluctuating ap- 
proximation do not affect the center line of the ribbon and 
hence do not decrease the mean square diameter. 



3 Results and Discussion 

To study the effect of a ring geometry on a fluctuating 
polymer ribbon with minor bending stiffness ai, relative 
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Fig. 2. Contour plot of the descent of the mean square diame- 
ter — d(D 2 } /d(Rc/ai) versus a — 01/02 and r = ai/03, where 
the contour lines decrease from n/2 at (00, 1) to 2/-7T at (0,0) 
in ten equal steps. The additional white contour marks where 
in-plane and transversal fluctuations are equal in decreasing 
the diameter. 



bending stiffness a = 01/02 and relative twist stiffness 
t = ai /a 3 we choose the mean square diameter as param- 
eter for the ensemble average since this parameter is ex- 
perimentally accessible and enables comparison with the 
mean squared end-to-end distance of an open polymer. In 
the limit of small fluctuations the analysis in the previous 
section enables us to evaluate the mean square diameter 
of a polymer ribbon constrained to a ring by inserting 
the modes of the Euler angles as given in equations (fT0| 
and (llip into the formula (fT5|) and calculating the sums, 
resulting in: 



(D ) = ARi 



ai 



tan(7r v / (l-a)(l-r)/2) 



2tt(t-1) ^(l-^Xl-^) 3 



(16) 



Analyzing this mean value the effects of the ring geome- 
try are unveiled, namely coupling of twisting and bend- 
ing and an effective stiffening of the polymer. To this end 
a plot of the diameter is desirable. Since that would re- 
quire a four dimensional space, we refrain from depict- 
ing the reciprocally proportionally decrease of the mean 
square diameter with bending stiffness a\ and show a 
contour plot of the descent of the mean square diameter 
—d(D 2 )/d(Rc/ai) versus a and r in Fig.[5J The diameter 
decreases by 2/n-Rl/ai at (r = 0, a = 0) up to 7r/2-i?j?/ai 
at (00, 1). 

Our analysis of the elastic free energy in the limit of 
small fluctuations shows a bend-bend coupling in addition 
to twist-bend coupling. The motion describing twisting Sip 
in equation ([3]) contributes to the transversal and the in- 
plane bending term in the elastic free energy. Vice versa 
the motion characterizing transversal bending 69 increases 
the twisting free energy term, hence bending and twisting 
are coupled. Furthermore, transversal bending motion in- 
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dicated by S9 in equation ([3]) decrease the in-plane bend- 
ing term in the elastic free energy, therefore both bending 
motions are coupled. In the case of the mean square diam- 
eter the coupling comes apparent in the complicated de- 
pendence on the bending and twisting stiffness in equation 
(fl6|) , in contrast to the uncoupled case, where the subtra- 
hend of the diameter is constituted of three independent 
term proportional to the inverse of the bending or twist- 
ing stiffness, respectively. Both transversal bending and 
twisting are strongly dependent on the in-plane bending 
stiffness, as all fluctuations are diminished for infinite in- 
plane bending stiffness a\. In the limit of a\/R c — > the 
motions become uncoupled as expected for the limit of a 
open polymer. The quality of the twist-bend coupling be- 
tween transversal bending and twisting can be read of the 
descent of the mean square diameter depicted in Fig. [51 
The coupling is strongest in the regime of small a and r, 
there a small increase in a is compensated for by a small 
decrease in r, yielding constant descent. The increased 
transversal fluctuations are reduced by a higher inabil- 
ity to twist. For larger a the coupling to r becomes less 
and less as an increase in a requires a larger and larger 
decrease in r for compensation. The twist-bend coupling 
vanishes for equal bending stiffnesses a = 1 or zero twist 
stiffness t = oo as expected. 

The effective stiffening due to the ring geometry can 
be analyzed by comparing the amount by which the in- 
plane bending modes and the transversal bending modes 
make the diameter decrease from its mechanical equilib- 
rium value. For an uncoupled open polymer ribbon both 
fluctuations are equal for equal bending stiffnesses. How- 
ever, if periodic boundary conditions are applied to gener- 
ate a polymer ring, in-plane fluctuations become strongly 
constrained while the transversal motion is only constrained 
via the bend-bend coupling. Discretizing a polymer ring 
to a polygon the restrictions on both bending motions 
due to the boundary conditions become obvious. While a 
single segment of an inextensible polygon can be moved 
transversal rotating by any angle up to 180°, a purely in- 
plane state is only achieved rotating by exactly 180°. All 
other in-plane bending moves would require neighboring 
segments to move as well, the motion is no longer truly 
local. In-plane bending is therefore strongly constrained, 
which decreases their fluctuation amplitude with respect 
to the transversal bending amplitude. Based on the ex- 
pression for the mean square diameter given in equation 
(fT5|) the origin of the subtrahends that make the diameter 
decrease of the fluctuating ribbon from D = 2R C and their 
amount are comparable. For r and a small the in-plane 
fluctuations account for the majority of the subtrahend, 
as transversal fluctuations are suppressed by large twist- 
ing stiffness and large transversal bending stiffness rela- 
tive to the in-plane bending stiffness. This changes if one 
of the relative stiffnesses increases, equal decrease of the 
diameter is reached along the white contour line in Fig. [5J 
Along this line ai is well larger than the in-plane bending 
stiffness ai, nevertheless transversal fluctuations extract 
more length from the diameter than in-plane modes. This 
puzzling behavior demonstrates the effective stiffening due 



to conformational constraints which mainly affects the in- 
plane motion. 

For reasons of illustration two limiting cases of the 
mean square diameter of a polymer ribbon are discussed 
in the following. First the limit of a double stranded bun- 
dle with a = is considered. Afterwards a polymer rib- 
bon with symmetric cross section, a = 1, is investigated. 
Two rigidly connected polymer strands suppress any fluc- 
tuations parallel to their interconnection, this behavior 
is recovered by our ribbon model in the limiting case of 
a-i — oo. Corresponding to a = 0, the mean square diam- 
eter is easily calculated to be: 



ai 



4 -\- tix 2 cot(7rVl - r) 



2tt(1 - t) 



tan(7T\/l - t/2) 
(1 - r) 3 / 2 



(17) 



Transversal fluctuations are still possible for any finite 
twist stiffness, they even increase rapidly from zero at 
t = as can be anticipated from the small distance of 
the contour lines for a = in Fig. [51 For r = oo the sub- 
trahend coincides with that of a = 1. Hence, zero twist 
stiffness 0,3 enables local twisting which washes out any 
difference between the bending stiffness ai and 02, even 
an infinite one. 

To verify our predictions based on the ribbon model 
by comparison with previous results and our simulations 
we study the special case of a symmetric cross section 
where both bending stiffnesses coincide ai = a 2 — l p , 
i.e. a = 1. This can be identified as an semiflexible poly- 
mer constrained to form a ring with persistence length l p 
and additional twist stiffness. The mean square diameter 
is computed to be: 



-1 = (2i? c ) 2 1 



1 2ttR c 



(18) 



This linear decay is up to one-hundredth identical to the 
numerical result obtained by Shimada and Yamakawa for 
the mean square radius of gyration [16j . which coincides 
with the mean square radius for semiflexible rings as the 
center of mass lies in the center of the rigid ring. Note 
further that the mean square diameter is independent of 
the twist stiffness, although no constraints where applied 
to the latter. In the limit of symmetric cross section the 59 
modes become independent of r although the Sip modes 
are still weakly coupled to bending, as can be seen by en- 
tering a = 1 into equations (fTTj) and (|T2|) . Since the mean 
square diameter depends on 6<f> and S9 modes only, the 
mean square diameter of a ring with symmetric cross sec- 
tion is independent of the twist stiffness. Comparing the 
mean square diameter in this limiting case with the end- 
to-end distance of a semiflexible polymer quantifies the 
effective stiffness resulting from the geometry of the ring. 
For a weakly fluctuating semiflexible polymer of length 
L and persistence length l p the end-to-end distance is 
(R 2 ) « L 2 (l — L/3lp). Examining the prefactors of the 
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L/l p 

Fig. 3. Comparison of simulation data (crosses) for the mean 
square diameter (D 2 )/R 2 versus L/l p to equation (|18[) (black 
line). For L/l p 3> 1 the mean square diameter decays with 
(L/l p )~ in accordance with previous results [16] (grey line). 
Error bars for the Monte Carlo data are approximately of the 
size of the symbols. 

linear term we note that a ring is effectively about five 
times stiffer than an unconstrained polymer. This effec- 
tive stiffening becomes also apparent when analyzing the 
shape of semiflexible polymer rings [25]. There, the stiff 
limit dominated by planar, elliptical shapes extends up 
to L/l p w 5, beyond which three dimensional, crumpled 
structures prevail. 

To asses the quality of the tight ring approximation 
a semiflexible ring with symmetric cross section and zero 
twist stiffness has been simulated using Metropolis Monte 
Carlo methods. The ring is described as a polygon com- 
posed of N tethers of fixed length a — (L/tt) sm(ir/N) and 
direction t|. The energy assigned to an individual config- 
uration is given by the elastic energy, F — YfcsT(Z p /L)x 

x X)i=i(l — t3*3 +1 )j imposing periodic boundary condi- 
tions , tjv+i = ti. New conformations are achieved by 
pivot moves [18], performing 10 6 Monte Carlo steps per 
segment. The results for the mean square diameter ver- 
sus flexibility L/l p presented in Fig. [3] suggest that our 
approximation is valid for L/l p up to seven, only at this 
flexibility the cross over to nearly flexible behavior as in- 
vestigated by Shimada and Yamakawa [TB] occurs. Note 
that there is no parameter to adjust. Since there is almost 
no cross-over region, a semiflexible polymer ring is entirely 
described by a stiff regime and an almost flexible regime, 
which are now both explained analytically. 
In general weakly bending approximations yield good pre- 
dictions for flexibilities up to one. However, our model for 
a weakly fluctuating ring is in agreement with simulations 
up to much higher flexibilities as a result of the effective 
stiffening due to the ring geometry. 



4 Conclusions 

In summary, we have investigated the effects of a ring 
geometry on a thermally fluctuating polymer ribbon. Pa- 



rameterizing the ribbon in terms of Euler angles and ap- 
proximating for small fluctuations an analytical expres- 
sion for the mean square diameter of the ring was derived. 
Analysis of this ensemble parameter and the elastic free 
energy of the fluctuating semiflexible polymer constrained 
to a ring reveiled coupling between in-plane and transver- 
sal bending motions as well as twisting motions. Further- 
more an effective stiffening of the ring mainly affecting 
the in-plane bending modes was found to result from the 
ring geometry and could be quantified as approximately 
five times the polymer's bending stiffness in the case of 
a symmetric cross section. Comparison with Monte Carlo 
data shows good quantitative agreement with our model 
up to high degrees of flexibility (L/l p « 7). Knowledge 
of an analytical result for the mean square diameter de- 
pending on two bending stiffness and twist stiffness en- 
ables new possibilities for the experimental determination 
of polymer and polymer bundle mechanics. As a thorough 
understanding of polymer properties is required for the 
unravel of cytoskeleton mechanics, we hope that our work 
will contribute to investigations in the broader field of cell 
rheology and bundle dynamics. 
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